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Molecular-dynamics simulations are presented for two correlation functions formed with the partial 
density fluctuations of binary hard-sphere mixtures in order to explore the effects of mixing on the 
evolution of glassy dynamics upon compressing the liquid into high-density states. Partial-density- 
fluctuation correlation functions for the two species are reported. Results for the alpha-relaxation 
process are quantified by parameters for the strength, the stretching, and the time scale, where 
the latter varies over almost four orders of magnitude upon compression. The parameters exhibit 
an appreciable dependence on the wave vector; and this dependence is different for the correlation 
function referring to the smaller and that for the larger species. These features are shown to be 
in semi-quantitative agreement with those calculated within the mode-coupling theory for ideal 
liquid-glass transitions. 
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I. INTRODUCTION 

If one compresses or cools a liquid, there appear slow 
dynamical processes which are referred to as structural 
relaxation. These processes are precursors of the liquid- 
glass transition. The study of these phenomena has been 
a very active field of research in recent years. Several new 
experimental techniques were introduced to measure the 
evolution of structural relaxation spectra within the GHz 
band. Molecular-dynamics simulation techniques have 
been improved considerably so that correlation functions 
of liquids in equilibrium can be obtained for time inter- 
vals covering five to six orders of magnitude. The wealth 
of information obtained on glassy dynamics is a challenge 
for the theory of amorphous condensed matter. However, 
consensus on the understanding of the slow dynamics in 
glass-forming liquids has not yet been achieved P, Q . 

Simple monoatomic liquids crystallize before structural 
relaxation dynamics is fully developed. Therefore, stud- 
ies of the glassy dynamics have to be performed on sim- 
ple molecular systems or suitable mixtures. Recently, 
for example, a four-component mixture was studied by 
neutron-scattering spectroscopy. This system transforms 
to a metallic glass at low temperatures, but it exhibits 
the same scenario for the evolution of structural relax- 
ation as known for molecules The first molecular- 
dynamics studies of structural relaxation in an equilib- 
rium liquid were performed for a binary mixture ofpar- 
ticles interacting by purely repulsive potentials |j, |^ lg . 
A binary Lennard- Jones system has been introduced Q , 
whose interaction potentials are similar to the ones pro- 
posed for the description of the glass-forming Ni-P mix- 
ture JJ . This system has been used extensively to analyze 
all facets of glassy dynamics in the equilibrium liquid and 
also for the quenched non-equilibrium system [^llCll|. In 
the mentioned previous studies, mixing was merely intro- 
duced as a means of suppressing crystallization. In the 
present paper, we analyze the influence of mixing on the 



structural relaxation. 

In order to identify the effect of mixing on the glassy 
dynamics, we have performed molecular-dynamics simu- 
lations for four binary hard-sphere mixtures differing in 
the size ratio of the constituents and in the composition. 
By increasing the total packing fraction up to 0.605, the 
evolution of structural relaxation was detected for a time 
interval up to five orders of magnitude. As reported ear- 
lier 1^ , two scenarios for mixing effects have been iden- 
tified. For a mixture with a small size disparity of the 
constituents, the increase of the mixing percentage of the 
small particles for a fixed total packing fraction leads to 
a slowing down of the long-time dynamics. In this case, 
mixing stabilizes the glass state. However, upon mixing 
particles with a large size disparity, the increase of the 
percentage of the small particles at fixed packing frac- 
tion speeds up the structural relaxation. In this case, 
mixing stabilizes the liquid. The present paper reports 
a detailed analysis of the long-time relaxation processes, 
traditionally referred to as alpha processes, for the two 
scenarios mentioned. 

Glassy dynamics and a liquid-glass transition can also 
be observed experimentally in colloidal suspensions. In 
particular, one can prepare glass-forming colloids where 
the interaction potential is a very good approximation to 
a hard-sphere repulsion |12| . Informative light-scattering 
studies of structural relaxation for such a hard-sphere 
suspension have been reported ■ To suppress crystal- 
lization, a narrow distribution of particle sizes was cho- 
sen. Strictly, such a system is a multi-component mix- 
ture. But ignoring the small polydispersity, it can be 
viewed as a one-component system. In the same sense, 
one can consider a colloid studied by Henderson et al. 0| 
as an approximation for a binary hard-sphere mixture 
with a size ratio 0.8 between the two groups of particles, 
and a mixture studied by Williams and van Megcn 15] as 
one with size ratio 0.6. For the first mixture, a dramatic 
effect of mixing on the nucleation ratio was observed, but 
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no effect on the glassy dynamics has been reported. For 
the second mixture, it was shown that the time scale for 
the alpha relaxation decreased upon mixing. Our sim- 
ulation results suggest that the cited experiments 
do not deal with colloid-specific features. Rather, they 
exemplify the two scenarios for mixing effects on struc- 
tural relaxation. Therefore, the present paper provides 
a detailed list of quantitative predictions for correlation 
functions of glassy colloids, which can be measured by 
photon-correlation spectroscopy. 

The mode-coupling theory for ideal liquid-glass transi- 
tions provides a physical explanation for the evolution of 
structural relaxation in simple systems and allows for the 
first-principle evaluation of the density-correlation func- 
tions [16^. The results of this theory for the hard-sphere 
system have been used for a detailed analysis of the light- 
scattering data obtained for hard-sphere colloids with a 
small polydispersity ^^l- The theory has been extended 
recently to a discussion of binary hard-sphere mixtures. 
In particular, the above mentioned two mixing scenar- 
ios had been obtained 0|. It was possible to describe 
a major part of the scattering data for a mixture 
quantitatively by the theoretical results • These find- 
ings provide a motivation to use our simulation results 
also for a detailed quantitative test of the mode-coupling 
theory for the alpha-relaxation process. 

The paper is organized as follows. In Sec.^ the simu- 
lation details are described, representative results for the 
two mixtures considered are exhibited, and the mode- 
coupling-theory formulas for the alpha-relaxation process 
are listed. Then, in Sec. IIIII it will be explained how 
the alpha-relaxation processes are parameterized. The 
results for the parameters are compared with the corre- 
sponding ones obtained from the mode-coupling-theory 
findings. Section Hvl presents a comparison of the alpha- 
relaxation master functions for the density-fluctuation 
correlation functions of the simulation data with the cor- 
responding theoretical results. In Sec.0 the findings are 
summarized. 



II. BASIC CONCEPTS AND RESULTS 

A. Specification of the systems 

Let N and g denote the total number of particles 
and the total number density, respectively, of the bi- 
nary hard-sphere mixture (HSM) to be studied. Further 
numbers specifying the system are the particle diameters 
da, the particle masses nia, their thermal velocities 
the partial number densities Qa and number concentra- 
tions Xa = Qa/ Q ^ Na/N, as wcll as the partial packing 
fractions (pa = {n /Q)Qa<fia- Here, a — A and B labels 
the big and small particles, respectively. In the present 
work, the size ratio 5 = d^/dp^, the total packing fraction 
Lp — if A + fB, and the relative packing fraction of the 
smaller species, x = Lp^/ip, will be used as convenient 
control parameters to characterize the thermodynamic 



state. 

Density fluctuations of species a for wave vector q are 
constructed from the positions r^"\ k = \,...,Na., of 

particles of type a: Qa{q) = '^k^^^[^^k'^]l ■ The 
partial structure factors Safsiq) — {Qa{<l)* Qf3{(t)) provide 
the simplest statistical information on the equilibrium 
distribution of particles, i.e., on the structure. Here, (•) 
denotes canonical averaging. The structure factors de- 
pend on the wave vector only via g = They can be 
written as Sap{q) = XaSap + QXahaf3{q)xi3, with h^piq) 
denoting the Fourier transform of the pair-correlation 
function. The latter can be expressed in terms of the di- 
rect correlation functions Ca/s^q) via the Ornstein-Zernike 
equation: ha/siq) = Capiq) + QY.^'^a^i.l)^iKp(l) ES- 
The Sap{q), haf}{q), and Ca[3{q) are elements of real 
symmetric two-by-two matrices. The discussions will 
be restricted to such stable and metastable states where 
Sapiq) and Cap{q) are smooth functions of q and of the 
control parameters p, S, and x. 

The main quantities of interest in this paper are the 
density correlators ^ai3iq,t) = {Qa{q,t)* Qisiq)) ■ These 
correlation functions provide the simplest statistical char- 
acterization of the structural dynamics. They are real 
even functions of time t, and they form the elements of 
a symmetric two-by-two matrix. In principle, the corre- 
lators can be measured as intermediate coherent scatter- 
ing functions by neutron-scattering experiments for con- 
ventional liquids, or by photon-correlation spectroscopy 
for colloidal suspensions. A short-time expansion yields 

+ 0{t^) |20|. 

Within the regime of normal-liquid states, say, Lp < 0.4, 
the short-time dynamics varies on a 30% level upon 
changes of the control parameters. There are no struc- 
tural relaxation phenomena apparent in the transient dy- 
namics. 

Two mixtures shall be considered in the following. A 
system with 6 = 0.60 and x = 0.20, referred to as the 
6 = 0.60-system, is representative for a mixture with 
large size disparity. In this case, 54% of the particles 
are of species B. A system with S — 0.83 and x — 0.37, 
referred to as the S — 0.83-system, contains 50% of small 
particles. It is representative for a mixture with small 
size disparity. These mixtures have been used before, 
together with systems of smaller percentages x, in order 
to demonstrate the evolution of mixing anomalies with 
variations of x • Since we are not interested in details 
of the short-time dynamics, the masses of the particles 
are chosen equal, i.e. va — vb- The units of length and 
time are chosen such that c^a = 1 and va — vb — With 
these units, the natural time scale for the microscopic 
motion is <niic = /^A = 1 • 



B. Results from molecular dynamic simulations 

We perform standard molecular dynamics simulations 
for binary mixtures of = 1237 and N ~ 700 hard- 
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sphere particles with size ratios S = 0.60 and 5 = 0.83, re- 
spectively. The algorithm follows the usual event-driven 
scheme for the simulation of hard-sphere particles [2ll |. 
where the trajectory of the system is propagated from one 
collision to the next one. To generate dense enough ini- 
tial configurations without particle overlaps, we applied 
the same procedure as described earlier |23 |: Starting 
from a random distribution of points, particles were sep- 
arated growing their diameters in successive steps until 
the desired size was reached. From the initial configura- 
tion, each simulation proceeds by an equilibration run, 
followed by a production run during which positions and 
velocities are saved for subsequent analysis. In all cases, 
the equilibration time was larger than the time it takes 
for the particles' average displacement to reach one diam- 
eter of the large species, cIa- Up to four independent runs 
per state point have been performed to reduce statistical 
errors. Density correlation functions $^^(5,^) and static 
structure factors Sap{q) have been calculated by averag- 
ing over the independent runs and over 300 different wave 
vectors q of the same modulus q. The longest simulation 
run requested about three weeks, i.e. the largest density 
studied took about three months of CPU time on a fast 
AMD Athlon processor to be completed. 

We checked that no crystallization occurred during 
the production runs by monitoring the time evolution 
of the pressure of the system, and by visual inspec- 
tion of the configurations. We also evaluated the wave- 
vector resolved structure factor without angular averag- 
ing to make sure that no crystalline peaks have devel- 
oped. Other mixture compositions than the ones pre- 
sented below have been tried; for 6 = 0.83 and x — 0.276 
as well as for S — 0.60 and x = 0.10, it was also possi- 
ble to study the glassy dynamics in the liquid phase , 
despite a stronger tendency to crystallization. A system 
with 5 = 0.60 and x = 0.05 did not stay in the homoge- 
neous liquid phase long enough for a study of structural 
relaxation. 

Figures n a-nd 121 exhibit a typical set of structure fac- 
tors and the corresponding pair distribution functions. 
The results refer to the S = 0.60 system, and the lines 
are calculated using the Percus-Yevick theory for HSM 
[23, 13 . Obviously, this approximation theory accounts 
for the data rather well, even though the packing fraction 
(fi — 0.60 considered is rather large. But there are small 
systematic deviations of the kind known from the discus- 
sion of hard-sphere mixtures at smaller packing fractions 
[2^. For example, the theory overestimates the height 
of the first and second peak of S'aa(9) by about 10%. 
The contact values for the radial distribution function 
are 11.6, 6.87, and 8.76 for the AA, BB and AB func- 
tion, respectively, while the Percus-Yevick theory yields 
8.88, 6.32, and 7.28, respectively. It will be discussed be- 
low that these discrepancies have to be acknowledged if 
one intends to consider the results of the mode-coupling 
theory (MCT) quantitatively. Even though the results of 
Percus-Yevick theory are well-known, a side remark on 
its qualitative features might be in order. Increasing the 
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FIG. 1: Partial structure factors Sapiq), a,l3 = A,B, of a 
binary hard sphere mixture (HSM) for diameter ratio of the 
particles 5 = de/dA = 0.60, a relative packing ratio x of 
the smaller species of 20%, and total packing fraction ip — 
0.60. Circles are the molecular-dynamics simulation results, 
while full lines are the results calculated within Percus-Yevick 
theory. Here and in the following figures, the diameter of the 
larger spheres is chosen as the unit of length, dA = 1. 



size disparity, i.e. decreasing S below unity, the height of 
the first diffraction peak in SAA{q) decreases. Simultane- 
ously, the wing of the peak at qdA ~ 9 increases. Within 
MCT, the first trend stabilizes the liquid state, while the 
second trend stabilizes the glass. The first trend domi- 
nates at the glass transition for small S, the second one 
that for larger S ^^l- For the radial distribution func- 
tions, these trends correspond to a reduction of the con- 
tact values and a decrease of the averaged radius of the 
first neighbor shell with decreasing 5. Moreover, at larger 
size disparities, a splitting of the first-neighbor-shell peak 
into a double peak is observed, as is evident in Fig. 13 

Figures |31 and 0] exhibit representative examples for 
the evolution of the glassy dynamics upon compressing 
the S = 0.60 mixture. The wave vectors qdA = 5 and 
qdA — 10 have been chosen since they exhibit the char- 
acteristic differences in mixing effects for small and large 
wave vectors that have been discussed before |0, [l8l| . 
Note that the values of the partial correlation functions 
can be quite different; in particular the small values of 
the AB correlator at qdA = 10 are responsible for the 
worse signal to noise ratio observed there. 

The correlators for = 0.40 are close to exponentials 
whose characteristic decay time is near the natural time 




FIG. 2: The pair distribution functions gafiir) corresponding 
to the results shown in Fig. 



scale imic — 1 for the microscopic dynamics. They are 
typical for normal-liquid behavior which can be described 
on a 30% accuracy level by Enskog's theory for dense 
gases 26j. If the packing fraction increases to = 0.45, 
the time scale for the decay of the correlators increases by 
about 40%, so that for t « 5 all correlators have decayed 
to below the 5% level of their initial values. Increasing (/? 
above 0.45, a new relaxation pattern evolves for the dy- 
namics outside the transient regime, say for < > 5. The 
^ai3 {q, i)-versus-log t diagrams exhibit a two-step relax- 
ation scenario that has repeatedly been observed before 
in simulation studies and experiments. First, the correla- 
tors decrease towards some plateau. The curves become 
flatter and the plateau lengths increase if tp increases. 
Then, the correlators decrease from the plateau to zero. 
The dynamics for t > imic is called structural relaxation. 
Our simulations document this process, which is charac- 
teristic for glass-forming liquids, for a time interval ex- 
tending over nearly five orders of magnitude. 

The second step of the structural relaxation, i.e., the 
decay below the mentioned plateau is conventionally re- 
ferred to as the alpha process. The figures demonstrate 
that the time scale for the alpha process increases the 
faster with increasing ip the larger the packing fraction. 
The decay cannot be described by an exponential func- 
tion, rather it is stretched over wider time intervals. Ob- 
viously, the whole structural relaxation pattern, in par- 
ticular the alpha process, shows a subtle dependence on 
the wave vector q. It is the goal of this paper to charac- 
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FIG. 3: Molecular-dynamics-simulation results for the den- 
sity correlators "l>a^(g, f), a,/3 = A,B, for the S — 0.60 mix- 
ture for wave vector qd-A = 5. The dotted and dashed lines 
refer to packing ratios </? = 0.40 and 0.45, respectively. The 
full lines are correlators for = 0.530, 0.550, 0.570, 0.580, 
0.590, 0.595, 0.600, and 0.605 (from left to right). Here and 
in the following figures, the unit of time is chosen such that 
the thermal velocity of the particles is unity. 

terize the alpha process, in particular its q dependence, 
quantitatively. 



C. Some mode-coupling-theory equations 

Within MCT, the concept of a plateau and of an alpha- 
relaxation process can be defined precisely in the sense 
of asymptotic laws describing the dynamics near an ide- 
alized liquid-to-glass transition. These laws provide a 
motivation for the parameterization of the data. In ad- 
dition, our data shall be used to test quantitatively the 
results of the theory. In this section, the required formu- 
las are compiled. 

Let us introduce an obvious matrix notation to get 
the following equations in a transparent form. S{q), 
^{q,t), etc., shall denote two-by-two matrices with el- 
ements Saiffiq), $a/9('7, i)j etc. The Zwanzig-Mori for- 
malism 20] can be used to derive the exact equation of 
motion for the density correlators, 

I{q)d^^{q,t) + S-\q)^{q,t)+ f M{q,t-t')dt''P{q,t') dt' 
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FIG. 4: The analogous results as shown in Fig. |3 but for 
wave number qdA = 10. 



In the sum, the abbreviation p = q — k is used. The 
vertices Vaa'a"{Q, k,p) are expressed in terms of the di- 
rect correlation functions Cap{k) and a triple average 

{Qa{q)* Qfi{k)Q-^{p)) ■ Simplifying the latter by the con- 
volution approximation, one gets 

"(Qt^jP) — {qk/q)Caa'{k)Saa" +{<}P/q)Caa"{p)Saa' ■ 

(Id) 

Note that the mode-coupling functional is determined by 
the equilibrium structure alone. Specifying the struc- 
ture factors and a model for 7Vf'°^(q,t), the preceding 
equations are closed. To proceed towards a numerical 
solution, one introduces a grid of equally spaced wave 
numbers extending up to a cutoff q*. In this paper, we 
use 140 wave numbers and qd^ — 56.8 for the calculations 
based on the simulated structure factors. Some results 
based on the Pcrcus-Yevick approximation shall also be 
shown, and for those we used 200 wave numbers up to 
qd\ = 79.8. We refer to Ref. and the papers quoted 
there for further details. 

The MCT equations exhibit bifurcations for the long- 
time limits of the solutions. For packing fractions ip be- 



Here, I{q) is a matrix of inertia parameters, Iap{q) = 
^1 {.<fv1^Xa)5af}- The kernels M^piq^t) are fluctuating- 
force correlation functions, and they reflect the com- 
plicated many-body interaction effects. The equation 
of motion has to be solved with the initial conditions 
*(g,i = 0) = S{q) and dt^{q,t = 0) = 0. The essen- 
tial step in the theory is the application of Kawasaki's 
factorization approximation in oder to identify the ker- 
nel contribution M^'~^'^ {q,t), which expresses the cou- 
pling of the forces to the density fluctuations. The re- 
mainder of the kernel, M"^(g, t), is assumed to describe 
the normal-liquid-state dynamics. It is anticipated to 
vary regularly in the control parameters and to decay 
on the scale tmic for the transient motion. One gets 
M{q,t) = M'"'^{q,t) -f M^^'^(g,i), where the mode- 
coupling kernel is given by the mode-coupling functional 



M^^^(g,i)=^[*W](a), (lb) 



which reads 



(Ic) 

I 

low some critical value Lp'^ = Lp'^{5,x), one gets ^{q,t — s- 
oo) = 0. In this parameter regime, the solutions describe 
ergodic liquid dynamics. For ip > ip"^, the long-time limits 
are non-degenerate symmetric positive definite matrices 
^{ijt oo) = F{q). For these states, the solutions de- 
scribe amorphous solids, i.e., ideal glasses. The long-time 
limits obey the implicit equations 

F{q) = Siq)-[Siq)-'+:F[F]{q)]-\ (2) 

These equations are defined by the equilibrium structure 
alone; neither the inertia matrix I(q) nor the regular 
memory kernel M'"^^{q,t) enter. The above equation for 
ffq ) can be solved by a standard iteration procedure 

Let F'^iq) denote the non-degenerate positive definite 
matrix of long-time limits at the transition point ip — 
if'^. For reasons of continuity, $(<7,t) has to be close 
to F'^{q) for a large time interval if \lp — (p"^] is small. 
The correlators are the closer to F'^{q) the smaller \(p — 
ip'^l, and the time interval of this close approach extends 
simultaneously. Thereby, the evolution of the plateaus, 
which were discussed above in connection with Figs. |3| 
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andU is explained by MCT, and the F'^iq) are the MCT 
expressions for the plateaus. 

The decay of the correlators from the plateaus to zero 
for small negative ip — ip'^ shall be characterized by some 
time scale T{(p). Obviously, Ihn^^ipa T{(p) = oo. Let us 
consider the dynamics of the relaxation from the plateau 
on the time scale r by writing t ~ tr with t fixed but 
positive. There holds 

lim^*(g,?T) = *(q,f), (3) 
where f) obeys the equation 

<ifiq,i) = Siq)Miq,i)S{q)-Siq)^ f M{q,i-i')^{q,i' 

at Jq 

(4a) 

to be solved with the initial condition ^{q,t = Q) = 
F'^{q)- Here, M{q, i) is determined by the mode-coupling 
fimctional for the critical point, 

M{q,i)^J^-[^{i)]{q). (4b) 

The numerical solution of the equation for "^{q, t) is done 
similar to that for the full equations of motion. 

Equation ^ implies the following conclusion. Given 
some and some error margin. 



from the simulation work. For the two mixtures we have 
calculated 

<^^(^^ = 0.548, 6 = 0.44 (5 = 0.60, x = 0.20); (7a) 
<^^(.^ = 0.545, 6 = 0.43 (J = 0.83, x = 0.37). (7b) 

From the simulation data, the critical packing fractions 
Lp'^ for the liquid-glass transitions of the two mixtures 
have been determined from the alpha-relaxation times 
of the density auto-correlation functions, see Sec. IIII El 
yielding values of 0.606 and 0.586, respectively. The er- 
rors of 10% and 7%, respectively, exhibited by the values 
noted in Eqs. H7a|l and (|7b|) . indicate the uncertainty one 
should expect for MCT results. It is worth to stress that 
^pf one bases MCT on Percus-Yevick approximation for 
the structure factors, one gets as critical values for the 
two mixtures 0.520 and 0.515, respectively. Hence, the 
use of correct instead of approximated structure factor 
input to the theory improves the results for the critical 
points. It is remarkable that the modest errors of the 
Percus-Yevick theory, which are exhibited in Fig.^ lead 
to noteworthy changes in the MCT results for the critical 
points. 

III. PARAMETERIZATION OF THE 
ALPHA-RELAXATION PROCESSES 



*(g,i) =*(g,t/T) 



(5) 



A. Evolution of the alpha process 



is valid within the margin for t = t/r > provided 
^—^f is small enough. This is the superposition principle 
for the MCT alpha process. It describes the correlators 
in terms of a (/^-independent master function 'i'(g, t), and 
attributes the strong (/j-dependence to that of the scale r. 
Presenting the correlators as functions of f, the interval 
for t where they coincide expands to lower values of t if 
Yp — ip'^\ decreases. The master functions ^{q,t) depend 
only on the equilibrium structure. Neither the inertia 
parameters in I{q) nor the regular kernel M'°^(q, t) have 
any influence on \I'(q, t). These quantities enter the time 
scale T only. 

There are complicated but straight-forward formulas 
to evaluate from JF^ the so-called von Schweidler expo- 
nent fe, < 5 < 1, a critical amplitude H^^\q)^ which 
is a positive deflnite matrix, and a correction amplitude 
H^'^\q) (22,]. These quantities determine the von Schwei- 
dler expansion of the master functions. 



<b{q,t) = F^{q)-H'^^\qr + H'-'>{q)t 



r(2)/ 



(6) 



Here, terms of order have been dropped. Thereby 
an equation is obtained for the beginning of the alpha 
process. 

The MCT equations have been studied before for bi- 
nary HSM using the Percus-Yevick approximation for the 
structure factors the present paper, results will 

be presented using the structure factors S{q) obtained 



The evolution of the alpha-relaxation scaling law is ex- 
amined in Fig. 13 The upper panel is a typical example 
for the majority of correlators obtained in our simula- 
tions. For every packing fraction Lp, some time scale 
T = t{lp) can be defined so that the long-time parts of 
the correlators coincide if these are considered as func- 
tions of t = t/r. This coinciding part provides the master 
functions "^apiq^i) for the alpha process of the fluctua- 
tion considered. Upon increasing Lp, the t interval where 
Eq. © holds, expands to lower values of the rescaled 
time t. Thus, the observed scenario conflrms the MCT 
prediction. However, our data also exhibit violations of 
the above described scenario, which cannot be under- 
stood in the framework of MCT. These occur only for the 
BB correlators of the 5 — 0.60 mixture for wave vectors 
around the structure-factor peak position, qdx ~ 7; and 
this only for the two largest packing fractions examined, 
tp — 0.600 and 0.605. The lower panel of Fig. [Slshows a 
representative example. The following analysis of "^{q, t) 
shall therefore be based on those data sets that do not 
exhibit the described phenomenon; i.e. for ip < 0.60. 

The stretched exponential. 



<j>^iq,i) = ^„(g)exp Ui/Uq))M,) 



(8) 



is an often used empirical function for the description of 
alpha processes. It was introduced by Kohlrausch for the 
description of dielectric relaxation data. The description 
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FIG. 5: Density auto-correlation functions 'J'aa(5, and 
$BB(g, t) of the 5= ds/dA = 0.60 system for wave number 
qdA = 7.5, presented as functions of the rescaled time t — t/r. 
The upper panel exhibits simulation results for packing frac- 
tions ip = 0.605, 0.60, 0.595, 0.59, 0.58, 0.57, and 0.53; the 
lower panel for = 0.605, 0.60, 0.595, 0.59, and 0.58 (from 
left to right). For each (p, the scaling time r is chosen such 
that the curves coincide for large times. 



of the alpha-process master function by an amplitude — 
also called plateau value — Aa (q) , a time scale fa (q) , and 
a Kohlrausch exponent (3a{q) shall be used here as well. 
The dashed lines in Fig. |21 exhibit representative exam- 
ples for an analysis of the normalized auto-correlation 
functions 4>\{c[,t) — ^xp^{q,t)/ SxA{q)- The figures con- 
tain rescaled data for ip = 0.595 and 0.590, in oder to 
identify a major part for the interval of rescaled times t, 
for which the superposition principle is valid. Similarly, 
Fig. [71 exhibits examples for a fit of the stretched expo- 
nentials to the numerical solutions of the MCT equations. 
The q vectors are as close to the ones of Fig.^las permit- 
ted by the use of discrete wave vector grids. Note that 
the choice of the overall time scale is irrelevant for the 
discussion of the master functions 

Figures and confirm an observation often made 
earlier: the fits by Eq. ^ provide a good description 
of a major part of the alpha process. However, there 
are also systematic deviations between the fit function 
i) and the master functions "^aa{q, i)- This holds in 
a similar manner for the fits to the data and to the MCT 
results. The fit contains unavoidable systematic errors, 
because the fit parameters depend somewhat on the time 
interval chosen for the fit optimization. In our analysis. 
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FIG. 6: Normalized large-particle density auto-correlation 
functions (j>A{q,t) — ^AA{q,t) / SAA{q) obtained from the 
simulations for three wave numbers q and packing fraction 
ip = 0.595 (left full line) for the 6 = 0.60 mixture. The right 
full lines are the simulation results for ip = 0.590, rescaled 
on the first curve for large times t. The dashed and dash- 
dotted lines are fits to the alpha process by the Kohlrausch 
function, Eq. (|HJ, and the von Schweidler expansion, Eq. JSJ, 
with b = 0.44, respectively. The fit interval for the latter is 
indicated by the horizontal dotted arrow. 



the fit was done so that the large-t part is described best. 
Thereby, the errors of the fit appear solely for the small-i 
part of the master functions. 

Equation ^ suggests another fit formula, which is 
valid for the small-t part of the master functions. But 
the small-i part can be identified only to the limited ex- 
tend to which the scaling interval can be established. 
Consequently, also the fit using the von Schweidler se- 
ries contains unavoidable uncertainties. The dotted lines 
in Fig. exhibit representative examples for such fits. 
Fits could be achieved with von Schweidler exponents b 
chosen between 0.40 and 0.50. Therefore, the predicted 
exponents, Eq. 0, are confirmed within an uncertainty 
of ±0.05. All the results shown are obtained with the 
cited theoretical exponents. The remaining fit parame- 
ters shall be discussed below. 

Note that these results can depend somewhat on the 
time window chosen for the fit, which is 3.5 < t < 200 
for the fits discussed here. Similarly, the dotted lines 
in Fig. □ exhibit the results of Eq. ® for the MCT 
results. But here, the functions F''{q), H^^\q), and 
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FIG. 7: The analogous set of curves as in Fig. IHl but now 
obtained for the numerical solutions of the MCT equations. 
The packing fractions for the full curves are tp — 0.548 and 
0.545. The additional dash-dotted lines are the asymptotes 
given by the first two terms in Eq. (|SJ. 

H^'^\q), as well as the scale r are calculated from the 
MCT equations. In this sense, the dotted lines are not 
fit results. The discrepancies between the dotted and 
the full lines represent the ones between the full solu- 
tion of the MCT equations of motion and the specified 
second-order asymptotic description of the solution. It 
is reassuring that the discrepancies between the full and 
the dotted lines in Fig. exhibit similar trends as the 
ones shown in Fig. [7| A quantitative account of the dif- 
ferences between the results shown in the two figures is 
included in the discussions of the following sections. 



B. The plateau values 

The circles in Figs. |H1 and |51 exhibit the plateau val- 
ues of the two-step relaxation process, Aa{q)^ obtained 
by fitting Eq. Q to the simulation data for the normal- 
ized auto-correlation functions (fia{q,t). For qd^ < 6, 
the Aa{q) are very large; they almost reach their upper 
limit unity for q tending to zero. This is a typical mixing 
phenomenon. For one-component systems, A{q) « 0.4 is 
expected for small q [2^|33- The width qi/2 of the Aa{q)- 
versus-g curves, defined by Aa{qi/2) = 1/2, is about 7% 
larger for a = A in Fig. |51 than in Fig. O This means 
that the large particles are better localized in the mixture 




FIG. 8: Circles and squares represent the Kohlrausch am- 
plitudes A{q) determined by fitting Eq. ||HJ to the simula- 
tion data and MCT solutions, respectively, for the normalized 
auto-correlation functions ^aa{q,t)/ Saa{q) of the big parti- 
cles, a = A (top panel), and the small particles, a = B (bot- 
tom panel), of the 5 = 0.60 mixture. The full and dashed lines 
show the MCT plateau values faail) = F^dq) / Saa{q) calcu- 
lated with the simulation results for the structure factor and 
the Percus-Yevick approximation for the Sapiq), respectively. 




FIG. 9: Analogous results as in Fig. |H| but for the mixture 
with size ratio S = 0.83. 
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FIG. 10: The circles exhibit the plateau values fafj{q) for 
the 5 = 0.60 mixture, obtained by fitting the simulation 
data for the normahzed correlators 4>a/3(<?)/V Saa{q)Si3p{q) 
to the von Schweidler expansion, Eq. @, see text, a = A 
refers to the big particles, a = B to the small ones. The 
fit was optimized to the data for 95 = 0.59. The full lines 
are the MOT results for the normalized plateaus fa/iil) ~ 
^0/3(9) / \/Saaiq)SBi3{q)- FoT o. — fi , they agree with the lines 
shown in Fig. |H| The squares are fit results optimized for the 
data for 1^ = 0.595. 



with the larger size disparity. For the smaller B particles, 
the opposite trend is observed. Similarly, the plateaus for 
the A correlators exhibit some small peak near the posi- 
tion qdx = 7.5 of the structure-factor peak, while the B- 
correlator plateau values exhibit a pronounced minimum 
there. For the 5 = 0.60 mixture, Ax{q) exhibits a min- 
imum near qdp^ = 11, which is accompanied by a maxi- 
mum near qd^ = 12. Instead, the 5 = 0.83 system has a 
shoulder for Aa(<z) for 9 < g < 12. All these details are 
reproduced semi-quantitatively by the results obtained 
from the MCT values, which are shown as squares in the 
figures. There are only a few cases where the plateau 
values deduced from the data differ by up to 10% from 
the ones deduced from the MCT results: the B plateaus 
for the 5 = 0.60 system for qdx ~ 16 or the A plateaus 
for the (5 = 0.83 system for qdp^ ~ 17, for example. 

The symbols in Figs. E| and ^] show the re- 
sults for the plateau values fapiQ) obtained from fit- 
ting Eq. © to the simulation data, normalized by 
^0/3(9, t)/ \/ Saa{q)Spf3{q). The fit results obtained for 
the 5 — 0.60 mixture for ip — 0.590 and those for 
(/? = 0.595 differ by about 5%. This difference thus ap- 
pears as the inherent uncertainty of the data analysis. 
No such difference could be identified for the 5 = 0.83 
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FIG. 11: The analog of Fig.[Tni but for the & = 0.83 system. 
The fits are based on the data for Lp ~ 0.582. 



mixture. The full lines in the figures exhibit the MCT 
plateau values /;j_a((?) = F^p{q)/ ^ Saa{q)Sf3i3{q). For 
the diagonal functions, a/3 = AA and a/? — BB, these 
lines are identical to the full lines shown in Figs. |S1 and 
El Those lines should help to compare the plateau re- 
sults Aa{q) with the corresponding fit results for the 
faail)- Obviously, all qualitative features of the plateau 
fits based on the Kohlrausch function agree with the ones 
based on the von Schweidler series, both for the fits to the 
simulation data and for those to the MCT results. The 
MCT results for /aa(9) the S ~ 0.83 mixture are in 
perfect agreement with the simulation data, while /bb (q) 
is underestimated systematically by MCT. But the dif- 
ference is only about 5%, except for qdA ~ 14, where 
the discrepancy reaches about 10%. The deviations for 
/ab(9) are of similar size. For the system with larger 
size disparity, the discrepancies between data and MCT 
result is somewhat larger, but it is not seriously larger 
than the inherent uncertainty of the fits. 

Figure [7| exhibits also the leading term of the 
von Schweidler expansion. In general, accounting for 
the next-to-leading term of 0{P^) increases the range 
of validity of the von Schweidler expansion dramatically 
t28]. Indeed, a data analysis with a g-independent ex- 
ponent b is possible only if the 0{t^^) term is included 
|31| . However, if ip is not close enough to (p"^, it may hap- 
pen that the 0{t^^) terms cancel ©(i^**) contributions. 
In such case, the fit range may shrink upon inclusion 
of the 0{P^) terms [s^. This accident is demonstrated 
by Fig. \7\ for qdA = 10.2. Such phenomenon cannot be 
foreseen in an unbiased data analysis, which then, nec- 
essarily, must lead to errors in the fit amplitudes. This 
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FIG. 12: MCT plateau values f^i^q) for the one-component 
hard-sphere system, calculated with the simulated structure 
factor (solid line), with the Percus-Yevick approximation 
(dashed line), and with the Verlet-Weis expression for S{q) 
(dash-dotted line). The values for qd < 5 based on the simu- 
lated structure factor are unreliable and have been cut off in 
the figure. 



explains why the sign of the correction amplitude iden- 
tified in the lower panel of Fig. |H| differs from the one 
shown in the lower panel of Fig. [7| 

An obvious source of errors in the MCT results is due 
to using incorrect equilibrium-structure information in 
the mode-coupling functional. It was mentioned in con- 
nection with Eq. ([TJ that replacing the structure factors 
by their Percus-Yevick approximations increases remark- 
ably the difference between the MCT results for the crit- 
ical packing fractions and the results derived from the 
simulation data. The dashed lines in Figs. |S1 and |51 exem- 
plify the same phenomenon for the plateau values. These 
are the MCT results for the f^ail) based on the use of 
the Percus-Yevick structure factors. This approximation 
for the equilibrium structure leads to underestimations 
of the plateau values by more than 10%. Figures |H1 and 
1^1 demonstrate that MCT is so sensitive to the small de- 
viations between Sap{q) and the Percus-Yevick results 
which are exhibited in Fig. ^ that they lead to serious 
flaws in the quantitative predictions for the dynamics. 

The difference between the MCT results for the f^p(q) 
based on the Percus-Yevick approximation and those us- 
ing Safsiq) as obtained from numerical simulations is of 
similar size for both the S = 0.60 and the S ~ 0.83 
system. It is not a specialty of mixtures, as is demon- 
strated by Fig. El There, the comparison is repeated for 
the one-component hard-sphere system (HSS). The MCT 
result for the critical point based on the Percus-Yevick 
structure factor is tppy = 0.516, and the corresponding 
plateau values fpyil) shown as a dashed line. The 
Verlet-Weis modification of the Percus-Yevick structure 
factor |33j is an empirical improvement of S{q), in par- 
ticular for the contact values. Using this as input for 
the MCT, one gets as the critical point (pyw = 0.525. 
The corresponding /vw(9)' which are shown as the dash- 
dotted line, are systematically larger than the ones based 
on the Percus-Yevick results for S{q). We could obtain 



simulation data for the structure factor of the metastable 
HSS for packing fractions up to = 0.54. Beyond this 
packing fraction, crystallization was always taking place 
before particles could diffuse one nearest neighbor dis- 
tance, making it impossible to generate data meeting 
our equilibration criteria. Extrapolating the smoothly 
(/^-dependent S{q) wp to ip = 0.55, we calculated as the 
critical point <^hss ~ 0.546 and obtained the full line 
as the plateau values. It is remarkable that the differ- 
ences between f^iq) and fy^{q) for qdA > 10 are larger 
than those between fy^J^f{q) and fpyiq)- Note that we 
find, using the simulation results for the structure factors, 
(p'^{S — 0.83) < <Phss < 'p'^{S — 0.60), i.e., the system with 
small size disparity shows a change of upon mixing 
that is qualitatively different from the one seen for large 
size disparity. This MCT result is qualitatively the same 
as predicted originally on the basis of the Percus-Yevick 
approximation 18] and is confirmed by our simulations 



C. Von Schweidler-expansion amplitudes 

The amplitudes H^^{q) in Eq. ® are the most im- 
portant parameters quantifying the dynamics in a time 
interval where the correlators are close to their plateaus. 
The upper panel of Fig. lTSI exhibits a set of representative 
results. It compares the amplitudes ha\q) of the nor- 
malized master functions ^ aaiq,t) / Saaiq) obtained from 
the analysis of the simulation data for the S = 0.60 mix- 
ture with the corresponding quantities calculated within 
MCT. For the quantitative comparison, a scale factor S 
has to be adjusted since the arbitrariness of the time 
scale r implies an arbitrariness in the prefactor of 
ha\q) = Haa{q)/ Saa{q)- The /iq ''-versus-q curves ex- 
hibit a subtle structure. While h^^\q) has a maximum 
for qdA ~ 5.7 followed by a minimum for qd^ ~ 7.4, 
i.e. for a wave vector near the structure-factor peak po- 
sition, /ig ^ (q) increases monotonically to a maximum for 
qdA w 7.4. While h'-^^q) increases for qdA > 7.4 mono- 
tonically to a maximum for qdA ~ 11, ^3^(9) exhibits 
a sharp minimum for q slightly above 7.4 before it also 
reaches a maximum for qdA^n. For qdA> H, h'^\q) 
decreases monotonically, while h'^^ (q) has a minimum for 
qdA ~ 12.7 and then exhibits a broad maximum. These 
features are reproduced by the MCT results. The MCT 
results agree with the data on a 10% level, except for 
the B amplitudes for qdA near 10, where there are 20% 
discrepancies. 

(2) 

The amplitudes ha (q), which describe corrections 
to von Schweidler's law '^aa{q,i)/ Saa{q) — faail) ^ 
ha\q)P', exhibit a zero at some wave vector g*. For 
q < q^i the amplitudes are negative, and for q > q^ 
they are positive. These features and also the value 
(7b w 6 are reproduced by the MCT results. Notably, 
the MCT results for qdA < 12 still share some quali- 
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FIG. 13: Amplitudes /iL^'^'(g) = H2'a\q) / Sca{q) for the 
von Schweidler-law expansion of the normalized density auto- 
correlation functions aa{q,t) / Saaiq) for the 5 = 0.60 mix- 
ture. The dots and circles are the results of fitting Eq. @ 
to the simulation data for a = A (large) and a = B (small 
particles) , respectively. The full and dashed lines are the cor- 
responding results calculated from the MCT equations. Scale 
factors 5* = 15 and 5*^ are used for the comparisons of the 
h^a\q) and h^a\q), respectively, to account for the different 
time scales r used in the analysis of the simulation data and 
of the calculations. 

tative features with the results obtained from the fit to 
the simulation data, e.g. the sharp peak followed by a 
sharp minimum in /ig (g) at qdp^ fa 7, and the peak in 

(2) 

hj^ {q) at qdA ~ 11. Otherwise, one notices serious dis- 
crepancies between data and MCT results. For example, 
MCT predicts a particularly large range of validity of 
von Schweidler's law for the density fluctuations of the 
large particles for a wave vector qdA w 15. But the data 
analysis is done best for this case by using a correction 
amplitude near 0.2. The identified discrepancies signal- 
ize the limitations in the determination of a correction 
amplitude for an asymptotic law from data which cannot 
be chosen sufficiently close to the singularity. 



D. The stretching exponents 

The exponent /3 in Eq. (jS]) provides a convenient overall 
measure for the alpha-relaxation stretching. It quantifies 
in an averaged manner deviations of the alpha-relaxation 
process from a Debye law, (j>{t) cx exp(— i/f). The lat- 
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FIG. 14: The stretching exponents (3a(q) obtained by fit- 
ting the Kohlrausch function, Eq. JSJ, to the alpha-relaxation 
master functions for the density auto-correlation functions of 
the two mixtures with size ratios S of the particles 0.60 and 
0.83. The arrows in the two panels mark the von Schweidler 
exponents 6, Eqs. l(7|l. The dots and circles are obtained by 
fits to the simulation data for the big particles, q = A, and 
the small ones, a = B, respectively. The full and dashed lines 
are the corresponding results obtained by fits to the solutions 
of the MCT equations. 



ter is the universal result for the dynamics of a variable 
coupled to a white- noise field. Figure H"!] shows that for 
2 < g < 6, the exponent decreases considerably with 
increasing q. There is no difference between the fluctu- 
ations for the large and the small particles in this wave- 
vector interval; and the stretching is larger for the system 
with smaller size disparity. For the wave vectors near the 
structure-factor-peak position, qdA ~ 7, the stretching 
of the B-fluctuations of the S — 0.60 mixture is much 
bigger than the one of the A-fluctuations: /3b ~ 0.5 
versus (3a ~ 0.75. There is an indication of the same 
phenomenon for the S = 0.83 mixture. For larger wave 
vectors. Pa is somewhat larger than /3b for the mixture 
with large size disparity. For the 5 = 0.83 mixture, (3a 
equals (3b for qdA > 8 within the noise of the data. 

For qdA < 6, the MCT results overestimate (3a{q) by 
about 10%. The qdA ~ 7 anomaly and also the large- 
q variation for the S = 0.60 mixture are described well 
by the theory. For the system with small size disparity, 
MCT overestimates the qdA ~ 7 anomaly, and there is a 
slight trend to underestimate (3a{q) for large q. 

For large wave vectors, the MCT alpha-process-master 
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FIG. 15: Alpha- relaxation time scales Ta{q) obtained by fit- 
ting Eq. (|SJ to the data (symbols) and MCT results (lines). 
The coding is the same as used in Figs. 1131 and 1141 A mul- 
tiplication factor s is applied to the MCT results in order to 
account for the different time scales relative to that used for 
the data. 



functions approach the Kohlrausch law. In this asymp- 
totic regime, /3q,((j) = 6 [sj]. Figure [T^ illustrates that 
this theoretical result is consistent with the simulation 
data. 



E. The alpha-relaxation time scale 

The fit of Eq. ((SJ to the long-time parts of the correla- 
tors yields the time scales fa{q) for the alpha processes 
up to an overall scale r. The results for the fits allow for 
a comparison of the scales for the fluctuations of the large 
particles with those for the small ones. They also allow 
to discuss the q dependence of the relaxation times. In 
order to compare the scales from the data with those cal- 
culated from MCT, one has to flt an overall scale factor 
s. 

Figure lTBl exhibits the results for the two mixtures. For 
qdA < 6, TA(q) = tb((z), and both scales decrease with 
increasing wave vector q. These features are reproduced 
by MCT, but the theory overestimates the time scales se- 
riously. For wave vectors near the structure-factor-peak 
position, ta exhibits a pronounced maximum, and tb 
has a sharp minimum. The ratios of the scales is about 
5 and 3 for the systems with large and small size dis- 
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FIG. 16: Ratio of alpha-relaxation times T"(g) (upper panel) 
and product of alpha-relaxation times with the large-particle 
diffusivity D — Da (lower panel) for the 5 — 0.60 system as 
functions of the packing fraction tp (see text). The lines are 
guides to the eye. 



parity, respectively. This feature is reproduced well by 
MCT. For the 6 = 0.60 mixture, Tj^{q) exhibits a maxi- 
mum for qdA ~ 13, while fQ{q) has a shoulder there. For 
qdA > 15, the relaxation times decrease with increasing 
q. The time scales for the large particles are somewhat 
larger than those for the small particles. These features 
are reproduced qualitatively by MCT, but the theory un- 
derestimates the time scales by a factor 2 to 3 for wave 
vectors above the structure-factor-peak position. For the 
6 — 0.83 mixture, the scale TA^q) exhibits a shoulder 
for 10 < g < 13 in accord with MCT. The time fB(g) 
exhibits a minimum for qdA ~ 12, while MCT shows a 
kink there. Again, MCT underestimates the Ta{q) for 
qdA > 10. 

Let us consider the variation of the alpha-relaxation 
time scales as functions of the packing fraction if. To 
this end, we have determined a time scale r" (g, (p) for this 
process by arbitrarily choosing (f>cia{q,t = T"{q, ip)) = 0.1, 
for those values of q where the plateau values are still 
appreciably larger than 0.1. Equation Q formulates the 
scale coupling of the MCT results. While the time scales 
of two alpha processes, say ti and T2, diverge for vanish- 
ing distance parameter s = {(p'^ — Lp)jpf^ the ratio ti/t2 
is a smooth function of e. For example, let t\ and T2 refer 
to the alpha processes of the correlator (pAA{q,t) of the 
S = 0.60 system for qdA = 10 and qdA = 5, respectively. 
If e increases up to about 0.05, ti/t2 decreases linearly 
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with e by about 4%. The simulation results behave sim- 
ilarly, as is shown in the upper panel of Fig. ^| If the 
packing fraction increases from 0.50 to 0.605, the alpha- 
relaxation time scales of the density-fluctuation correla- 
tors increase by more than three orders of magnitude. 
Still, the shown three representative ratios for the scales 
vary by less than a factor two. Hence, the scale-coupling 
prediction is verified on a 10% level for the three ratios 
shown in the upper panel of Fig. ll6l bv the open symbols. 
These examples are representative for density-fluctuation 
scales with intermediate and large wave numbers. If one 
of the wave numbers decreases to small values, the vio- 
lation of the scale coupling becomes larger, as is demon- 
strated by the full symbols in Fig. El The diffusivity 
D is proportional to the inverse of the alpha-relaxation 
scale T]j of the mean-squared displacement. Hence, a 
T • D-veis\is-(p diagram demonstrates the coupling of the 
scales for the processes described by r and that for the 
diffusivity. The lower panel of Fig. El shows the results 
for the 6 = 0.60 system. Here, r refers to the scale for 
the density fluctuations of the A particles, and the D are 
the simulation data for the A-particle diffusivities [lll |. 
The scale coupling holds for e > 0.05. However, ap- 
proaching the transition point more closely, contrary to 
the MCT results, the scale for the diffusivity decouples 
from that for the density fluctuations. The diffusivity 
does not decrease with decreasing e as strongly as l/r. 
The results of our simulations for the S = 0.83 system 
behave similarly. The described decoupling is in quali- 
tative agreement with the behavior found earlier for the 
simulation results of a binary Lennard-Jones system ■f] 
and for a model for water [35j- 

MCT predicts a power-law divergence in the asymp- 
totic limit of vanishing e for the common scale t in 
Eq. T oc e~'^. The exponent 7 is determined by the 
von Schweidler exponent b The value b « 0.44 used 
throughout the preceding discussions implies 7 ~ 3.0. 
Figure El demonstrates this property for the MCT re- 
sults for Tx{q) for three representative wave numbers q 
in form of a rectification diagram. In agreement with 
typical results for the simple HSS the asymptotic 
description holds well for e up to about 0.05, and there 
appear deviations if differs from Lp'^ by more than 5%. 
Analogous rectification diagrams for the simulation data 
are shown in Fig. 1181 Linear extrapolations of the data to 
large Lp for three wave numbers q yield estimates for the 
critical packing fractions: p'^{5 — 0.60) = 0.606 ± 0.001 
and ^"{5 = 0.83) = 0.585 ± 0.001. The rectification 
curves for the diffusivity relaxation times for both types 
of particles are included in the figure as D^^^-YeYSvs-p 
plots. The data seem to follow the power-law predic- 
tions, but lead to slightly different estimates of p'^, also 
depending on the species a — A,B. For the 5 — 0.60 
system, we get (/j^^ « 0.609 and (/sfj^^ ~ 0.619, while for 
the 5 = 0.83 system, ip^,^ « 0.588 and « 0.590. The 
decoupling of the diffusivity scales from the ones for the 
density fluctuations mentioned above yields this overesti- 
mation of the critical packing fractions. Let us emphasize 




FIG. 17: T~^^^ with 7 = 3 as function of the packing fraction 
ip. The r are the alpha-relaxation times for the large-particle 
density-fluctuation correlators <j)AA{q,t) calculated from the 
MCT equations, defined by <j)AA{q,T) — 0.1. The wave num- 
bers q are 10.2 (diamonds), 7.4 (squares), and 5.0 (circles). 
The dotted lines are the asymptotic laws t~^^"' oc {p'^ — p) 
with p''{5 = 0.60) = 0.548 and p\5 = 0.83) = 0.545. 



that the described estimations of p'^ have been done with 
the bias of a given exponent 7 = 3.0. An unbiased three- 
parameter fit of the scale as a function of p by the formula 
T = Ciip''" — p)^'' suffers from correlations between the fit 
parameters p'^ and 7. Indeed, such fits to the diffusivities 
of the two species lead to differing exponents 7 in 
disagreement with MCT. A similar result has also been 
found in a simulation of a binary Lennard-Jones mixture 
0- Fig. 1 181 however shows that on the basis of our simu- 
lation data, one can, given the restricted range of validity 
of the asymptotic law, not distinguish between these dif- 
ferent 7 values. Note that the largest discrepancy, either 
in 7 or in p'^ , emerges for the B particles in the 5 = 0.60 
system. Unbiased fits lead to larger uncertainties for p'^, 
since a decrease of the fit parameter p'^ can partly be 
compensated by a decrease of the fit parameter 7. One 
could get the crossing points of the D^^'^ -versus-(p curves 
closer to that of the r~^/'''-versus-(y9 curves in Fig. 1181 if 
one would use some 7' < 7. Such formulation of the de- 
coupling phenomenon is suggested by Fig. 1161 since the 
increase in 13 • r for p increasing above 0.58 can be fitted 

by [p"^ — p^Y , x = 7' — 7 < 0. 

A comparison of Figs. 1171 andll8lleads to two questions, 
which we cannot answer. Why is the range of the dis- 
tances £, where the power-law asymptote describes the 
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FIG. 18: r"^''^ with 7 = 3 (left scale) as function of (f. 
The T are the alpha-relaxation times for the large-particle 
correlation functions at wave number qdA = 10 (diamonds), 
7 (squares), and 5 (circles), defined by (^AA(g, t) =0.1. The 
upper and lower triangles are D^^^ (right scale) for the dif- 
fusivities of the small and big particles, respectively. The 
dashed lines are linear interpolations of the data for large tp. 

alpha scales, so much larger for the simulation results 
than for the MCT solutions? Why do the deviations of 
T^^l^ for larger e from the small-e asymptote have a dif- 
ferent sign for the simulation results than for the MCT 
solutions? 




10"' 10° 10' 10' 10^ 



FIG. 19: Alpha-master functions 00,^(5, f) = 

4>ap{ci,t) / fafi(ci) for the 5 — 0.60 mixture {a = A refers 
to large particles, a = B to small ones). The circles and 
triangles represent the simulation data for ip — 0.595 and 
0.59, respectively, scaled as in Fig. |S] The solid lines are 
MCT results for the master function ^ afi[q,t) / F^p[q), and 
a solution for a distance e = {tp ~ ip'^)/p'^ ~ —0.015 from 
the critical point. Time has been rescaled for all curves in 
order to match the alpha-time scale of the simulation data 
at qdA — 5, a/3 — AA. The results refer to wave vectors 
qdA ~ 5, 7.5, 10, and 12.5 (from top to bottom), where the 
curves for different q have been shifted vertically by 0.3 for 
enhanced clarity. 



IV. THE ALPHA-PROCESS SHAPE 
FUNCTIONS 

The preceding parameterization of the correlation 
functions shows that there is no universal master func- 
tion for the alpha-relaxation processes. For example, 
the stretching of the density-fluctuation auto-correlation 
functions for the large particles generally differs from the 
one for the small particles, and it depends on the wave 
vector of the fluctuations. It is a challenge for a micro- 
scopic theory to describe the alpha-process shape func- 
tions for different probing variables. Figures [T^ and QUI 
present our simulation results for the two mixtures under 
consideration for four representative wave vectors in com- 
parison with the MCT curves. All results are rescaled by 
the plateau values and by the time scales discussed in 
Sec. mil The data are shown for two packing fractions in 
order to document the asymptotic regime of validity of 
the superposition principle. The MCT master function is 
complemented by a curve for some small ip'^ — (p, in order 



to indicate the effect of preasymptotic corrections to the 
alpha-process asymptote. 

Figure El shows that, typically, the decay of the cor- 
relators of the 6 = 0.60 system from 90% to 5% of their 
respective plateau values is stretched on a time interval 
of about three decades. With two exceptions, this de- 
cay is described well by the MCT master functions for 
the alpha process. The first exception is the AA correla- 
tor for qdx = 7.5 for rescaled times t « 4 x lO'^. Here, 
the data decay less rapidly as the exponential long-time 
tail obtained by MCT. The second exception concerns 
the BB correlator for qdA — 10 and qdA = 12.5. Here, 
the data fall on the master functions only for such long 
times, where the correlator is below 70% of the plateau 
value. But, the MCT results for {ip - if") /if" « -0.015, 
which are added in rescaled presentation as full lines in 
Fig. El exhibit the same phenomenon. The reason for 
this is the large size of the critical amplitude /13 ^ [q) for 
these wave vectors, cf. Fig. El They cause particularly 
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FIG. 20: As in Fig. ^ but for S = 0.83 and ip « 0.582 
(circles) and (p ~ 0.575 (triangles). The solid lines are the 
theoretical master curve and a curve with e « —0.01. 



large preasymptotic corrections to the alpha-scaling law 
near the plateau. Thus, the specified exceptions are not 
a defect of MCT, but a confirmation of a subtlety of that 
theory. 

The test of the alpha-process shape functions for 6 = 
0.83, Fig. 1201 exhibits a series of problems. Deviations of 
the data from exponential decay for the very long rescaled 
times t Ki 5 X 10^ occur for qdA = 5 and qdA = 7.5 for 
all three correlators. The AB correlator for qdA — 12.5 
shows more stretching than the MCT solutions. Further- 
more, the BB correlators for qdA ~ 10 and qdA — 12.5 
miss the plateau. However, the latter is an obvious mis- 
take of the data analysis, which could be eliminated by 
correcting the plateau value. We did not carry out this 
correction in order to emphasize that 5% errors in the 
determination of the plateau values are almost unavoid- 
able in an unbiased data analysis. The reduction of the 
scaling regime for the BB correlator for qdA = 10 and 
12.5 occurs for the S — 0.83 system as discussed above 
for the 6 = 0.60 system; and it can be explained in the 
same manner. 

The most severe problem exhibited by Fig. [501 is the 
following. Even for e ~ 0.02, in the simulation, several 
correlators stay within a 5% interval around the plateau 
for time intervals as large as 2.5 decades or more. This 
holds, e.g., for the AA correlator for qdA = 7.5, for the 
BB correlator for qdA = 5, and for the AB correlator 
for qdA = 7.5 and 10. The full lines show that this fea- 
ture cannot be explained by MCT. Even if the distance 



parameter e = {(p'^ — ip) jipf is as small as 10~^, the calcu- 
lated correlators cross their plateaus much steeper than 
exhibited by the simulation data. 



V. CONCLUSIONS 

Molecular-dynamics simulations have been presented 
for two dense binary mixtures of hard spheres. One mix- 
ture deals with a size ratio (5 = 0.60 for the two particle 
species and the other with 8 — 0.83. The first system is 
representative for the situation where mixing stabilizes 
the liquid state, and the other for the one where mixing 
stabilizes the glass JJJ. The data demonstrate the evo- 
lution of a two-step relaxation scenario with increasing 
packing fraction Figs. |2| and ^ which is similar to the 
one detected previously for other systems. In this paper, 
a comprehensive analysis of the second relaxation step, 
usually referred to as the alpha process, is presented. It 
deals with the decay of the correlators from some plateau 
to zero. The process was identified as that part of the cor- 
relators exhibiting the superposition principle predicted 
originally by mode-coupling theory (MCT), Fig.[Sl This 
pattern is exhibited by all our simulation data except for 
the correlators for the small particles of the 8 — 0.60 sys- 
tem for wave numbers near the structure-factor peak po- 
sition, and this for the two largest densities p — 0.600 and 
0.605, Fig.|5l This violation of the superposition princi- 
ple might indicate that MCT ignores relaxation processes 
which become important close to the liquid-glass transi- 
tion point in this mixture. A possibility that we can- 
not exclude, however, are precursors of crystallization or 
phase separation. 

The simulation data for the alpha process have been 
fitted by Kohlrausch functions and by the von Schweidler 
expansion, Eqs. © and ©, respectively. These fits pro- 
vide two estimates of the plateau values. Usually, these 
estimates agree on a 5% accuracy level, Fig. |H1 Analyz- 
ing the MCT results with Kohlrausch functions, one gets 
plateaus which agree with the correct values within 5%, 
Fig. [7| This means, that both mentioned fit procedures 
yield reliable estimations for the plateaus, within the in- 
dicated uncertainty level. These plateau values exhibit a 
remarkable structure as functions of the wave number q. 
There are qualitative differences in the structure of the 
plateau functions referring to the A particles and the B 
particles. And there are quantitative differences between 
the results for the two mixtures. All these results are 
described quantitatively by the MCT results except for 
some rare cases, where data and theory differ by up to 
10%, Figs.EHIIl 

MCT requires the structure factors as input for the 
equations of motion. The MCT results reported in this 
paper are based on the simulation results for the stud- 
ied systems. Fig. ^ Replacing these structure factors by 
their Percus-Yevick approximation results, MCT still re- 
produces all qualitative features of the mentioned plateau 
functions. However, one systematically underestimates 
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the data. The error of the MCT results caused by the 
specified erroneous structure information can be as large 
as 20%, Figs. IHl and El These findings do not deal with 
mixture-specific effects. They apply also for the simple 
hard-sphere system. Fig. ^| 

The stretching of the alpha processes is parameterized 
by the Kohlrausch exponents [3a (q) , and the relative time 
scales are quantified by the scales Tq, (g) . These quantities 
vary with the wave vector and they depend on the index 
a for the species, a = A,B. MCT reproduces the Pail) 
reasonably, but there occur discrepancies up to 20%. The 
trends for the Ta{q) are reproduced by the theory, but 
there occur large quantitative errors. Figs. 1141 and 1 151 

Writing the von Schweidler expansion, Eq. JHJ, 
for the diagonal correlators in the form \['Qct(9,t) = 



one notices that h^a\q) defines a relative time scale 
fo.[q) = [h''it\q)/F^^{q)]-^l'' while specifies the 

shape. The amplitudes haHq) are reproduced reason- 
ably by MCT, but there occur errors up to 20%. MCT 
describes the trend of the q- and species-dependence of 

(2) 

the correction amplitudes ha'{q), but there are large 
discrepancies between data and theory. Fig. 1131 

The contradicting conclusions concerning the descrip- 
tion of the alpha process parameters arrived at in the pre- 
ceding two paragraphs indicate that the described prob- 
lems are ones of the fitting procedures. Indeed, Figs. 1191 
and 1201 show that MCT describes the alpha-process mas- 
ter functions well. 

Qualitative discrepancies between MCT and simula- 
tion data concern the ratio of the time scales T(g) for the 
alpha-relaxation processes for the density fluctuations of 
intermediate wave numbers q and the time scale td deter- 
mining the strong variation of the particle diffusivity D, 
D (X 1/t£). If the packing fraction of the 6 = 0.60 system 
increases from 0.51 to 0.605, logr(q) increases by 3.5, but 
logTD increases only by 3.0. The ratio T{q)/To cx T{q)D 
is practically constant if ip increases from 0.51 to 0.58, i.e. 
there is perfect scale coupling within this density inter- 
val. But increasing ip further, the ratio increases by up 
to a factor 4, Fig. ^] MCT overestimates the trend to 
particle localization near the glass transition point. This 
overestimation of td seems to be the reason why the cal- 
culated relaxation times T{q) exceed the data for small 
q, as is demonstrated in Fig. for qdA < 5. The longest 
time scale for g ^ is that for the collective diffusion 
process: T'^°^\q) oc l/((7^£''^°"(g)). The coupling of this 
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mode to the density fluctuations causes the divergence of 
T{q) for q ^ 0. One should expect that an underestima- 
tion of the tagged-particle diffusivity D implies the same 
mistake for the collective diffusivity D"^"^^. This explains, 
why the decoupling of T{q) and T(g') increases if q' de- 
creases to small values, as shown by the filled symbols in 
Fig.Hni 

The increase of the time scales r(g) with increas- 
ing packing fraction is described well by the asymptotic 
power law for the MCT results. Fig. ^| However, the 
mentioned scale decoupling implies that the extrapola- 
tion to zero of the (l/ri3)^/^-versus-(/3 graphs leads to an 
estimation of the critical packing fraction tp'^, which ex- 
ceeds the value obtained from the (l/T(g))^/'''-versus-(y9 
extrapolation by 0.6% Tl^. 

Finally, a feature of our simulation data should be em- 
phasized which concerns the time regime where the corre- 
lators cross their plateaus. It deals with times larger than 
the ones describing the short-time transient but preced- 
ing the regime of validity of the alpha-relaxation scaling 
law. Within MCT, this regime is described for large den- 
sities by the /3-relaxation scaling laws. In this respect, the 
MCT results for the hard-sphere mixtures behave as the 
ones for the simple hard-sphere system |0. Figure 1^ 
shows, however, that the correlators for the 6 = 0.83 sys- 
tem are close to the plateaus for time intervals exceeding 
the ones for corresponding MCT results by more than 
an order of magnitude. Hence, the /3-relaxation theory 
cannot account for the simulation data dealing with the 
plateau crossing. In that respect our simulation data for 
the hard-sphere mixture are also qualitatively different 
from the ones measured for quasi-bidisperse hard-sphere 
colloids 15] and from the simulation data for the binary 
Lennard- Jones mixture ISl. 
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